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Schulz and Shastry Reply:The previous Comment 
on our Letter Ref[l] (SS) solving exactly a class of two 
chain fermi systems, draws attention to earlier work in 
Ref[2] where many lattice chains are coupled together. 
The Comment makes several general remarks, emphasiz- 
ing the similarity of the resulting solutions, and expounds 
on their significance in the light of recent literature. We 
believe that the similar "look " of the solutions notwith- 
standing, we are dealing with different classes of models. 

The comparisions are facilitated by presenting the 
main points of our recent generalization Ref[3] of the re- 
sults of SS to the case of arbitrary number of chains. 
Firstly consider the continuum: let a — \,u represent 
the index of one of v chains of fermions Na in number, 
having coordinates Vi^a- The model solved is given by the 


hamiltonian H = ^ J2a i^f^. 

mentum Il^^a = Pi,a + J2j,l3 ia,f3 Va,/3{ri^a - 

object ^a,i3 is antisymmetric £,a,f3 — —£,f3,a 


where the canonical mo- 
Tj^fs), and the 
=sign(/3 


The gauge potential Va^pix) = 
Va,f3ix + L)is obtained from 
E„,Jx) = Ef- 


Va.fi{-X) 

EL4x) = 


-E„ 


a). 

= Vi3,a{x) = 

- Va.fiix) sat- 
isfying E^p{x) = ^/3,q(x) = -E^p{-x). A pseudo 
unitary transformation maps H into a "free model" : 
exp(iS')7Jexp(-iS') = H = \Y,a,iP'i,a^ where the gen- 
erator S = ^J2i,j,a,i3^a,f3Ea,i3{ri^a -"Tj^p). Imposing 
PBC's on the wave function^ = &cp{—iS)ip free{{ki^a\) , 
gives us conditions on the pseudo momenta fci^L = 
2-K rii^a + J2i3 ^a,i3 ^a,i}Np, whcrc the phase shift ^a^p = 

L Va^f3{x) dx. Specalizing to i^ = 2 recovers the prevous 
results in SS. 

Next consider a collection of v chains of 
fermions with a lattice Hamiltonian Hi = V — 

Sm,a[(6XP«'/'m,a)4„_„Cm+i,a -|- h.c], whcrC (/)m^a = 

J2i B^a,f3 J^a,f3i'm — Vfui^p. We perform a unitary 
transformation with a suitably generalized S — 

h^l,7n,a,P^a.,/3Ba,0il - m) ni^anm,(3, with Ba.p{m) = 

Bi3,a{m) = —Ba,[i{—in), where the burden of symmetry 
in exchanging both spatial coordinates and species in- 
dex is partitioned between -Ba,/3 and ^a.js- We find H[ = 
exp{iS)Hiexp{-iS) = V - J2a{T,lZ.\{cL,aCm+l^a + 
h.c] + [expiQaCL a'^i.a + /i.c] },where the twist an- 
gle 9q — J2p^a.i3^a,pNi3, thc phasc shift is Aq_^ — 
X]m=i ^a,fj{fn), and the number operator of fermions in 
chain (3 is iV^, which can be replaced by its eigenvalue 
in each sector, since it is conserved. Thus the model be- 
comes a twisted b.c. fermi chain with a given twist angle, 
and we can write down the equation governing the pseu- 
domomenta Lkj^a = Ztt/j^q, -|- 8q -I- {Bethe Phase Shift}. 
In obtaining this simple form of iJ/, we have chosen, as 
in SS the functions A,B {a, /3 in SS) to satisfy recursion 
relations -Bq,,/3(to+ 1) — i?Q,/3(7n) = Aa,fs{m) for 1 < to < 
i — 1. This together with the oddness of i?Q,/3 forces a 


spatial symmetry on A : Aq_^(— m) = Aa^/^im — 1). Once 
again, restricting to z^ = 2 recovers the results of SS. 

From the above discussion of the lattice model, note 
that the minimal m,odel of our class has Aa^fj = except 
for Aq_/3(0) = Aa.i3{—1) ^ 0, corresponding to the phase 
factor 4>m,a = Y.i3^a,pAa^p{Q)[nra,p + nm+i,!3\- The 
model of Rcf[2] corresponds to (j)m,a = const [r7,„i_Q,_i — 
nm,a+i\- While the chain index dependence can be 
obtained by specializing our general A, the issue con- 
cerns the dependence on the spatial index. In or- 
der to obtain the model of Ref[2], we must choose 
Aa,i3{m) —C0Tist5mfi5o,-i3,±i- This spatial dependence is 
inconsistent with the symmetry requirement noted in the 
previous paragraph, and hence we find that the m,odel of 
Ref[2] cannot he solved by our method. Morover, there 
are important dfferences in the boundary conditions in 
the chain direction, we have no need for any specific bc's 
whereas the Ref[2] seems to. Therefore the dependence 
of the twist angle Oq on the number of particles in the 
boundary chains are different from Ref[2]. 

To sumarize, the class of models introduced in SS con- 
tain considerable freedom in regard to the range and na- 
ture of the interchain coupling, subject however, to cer- 
tain spatial-index symmetries of the functions A, B. The 
models considered in Ref[2] do not satisfy our conditions, 
and we are unable to solve the same. Detailed examina- 
tion shows that the solutions differ : despite the similar 
look, we are dealing with different classes of morfeZsRef [4] . 
The details of the properties of the multi chain models 
reported here will be published in Ref[3]. BSS would 
like to thank Bell Laboratories Lucent Technologies for 
hospitality, where the work of Ref[3] was carried out. 
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